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The properties of three-jet events with total transverse energy greater than 320 GeV and individual jet energy 
greater than 20 GeV have been analyzed and compared to absolute predictions from a next-to-leading order 
(NLO) perturbative QCD calculation. These data, of integrated luminosity 86 pb _1 , were recorded by the CDF 
Experiment for pp collisions at y^s = 1.8 TeV. This study tests a model of higher order QCD processes that 
result in gluon emission and can be used to estimate the magnitude of the contribution of processes higher than 
NLO. The total cross section is measured to be 466 ± 3(stat.)lyQ 7 (syst.) pb. The differential cross section is 
furthermore measured for all kinematically accessible regions of the Dalitz plane, including those for which the 
theoretical prediction is unreliable. While the measured cross section is consistent with the theoretical prediction 
in magnitude, the two differ somewhat in shape in the Dalitz plane. 



In perturbative QCD, hard scattering of the 
constituent partons in the proton and antipro- 
ton results in events with large total trans- 
verse energy, J2 Et ■ Outgoing scattered partons 
hadronize and may be detected as hadronic jets. 
Three-jet events can be produced when a hard 
gluon is radiated from any of the initial, interme- 
diate, or final state partons in an event with two 
primary outgoing partons. 

We analyze here some properties of the cross 
section for three-jet event production in proton- 
antiproton collisions at the Fermilab Tevatron 
Collider at center-of-mass energy 1.8 TeV. The 
data, which were recorded by the Collider De- 
tector at Fermilab (CDF) [1], are compared with 
predictions by the first complete next-to-leading 
order (NLO) QCD generator, Trirad [2], for 
hadronic three-jet production at hadron colliders. 
We compare the measured and predicted absolute 
cross sections to test our understanding of the 
higher order QCD processes that result in gluon 
emission and to estimate the magnitude of the 
contribution of processes higher than NLO. 1 In 

1 While no quantitative estimate of the contribution of 
next-to-next-to-leading order processes to the cross sec- 
tion is available at this time, considerable progress has 
recently been made in calculating two loop 2^2 parton 
processes [3], important groundwork for the future. 



some kinematical regions, we provide a measure- 
ment of the cross section where the theoretical 
prediction is not reliable; this measurement may 
be a useful guide for theoretical calculations. We 
also compared the shapes of the measured and 
predicted cross sections when normalized, to ex- 
amine the sensitivity of the cross section to varia- 
tions in the value of the strong coupling, a s . The 
data sample corresponds to an integrated lumi- 
nosity of 86 pb -1 collected during the 1994-1995 
run (Run lb). 

A previous paper [4] examined a smaller 
dataset and was limited to a comparison with 
leading order theoretical calculations [5]. A sub- 
sequent analysis [6] compared a larger dataset to 
predictions from the HERWIG [7] parton shower 
Monte Carlo program and to the NJETS [8] 
leading order 2 — ► N parton-level prediction. 
The NLO calculation used here has the bene- 
fit of reduced renormalization scale dependence 
(and consequently lower systematic uncertainty) 
as well as a more reliable description of multijet 
production throughout phase space. This study 
expands upon the previous investigations by com- 
paring the data to absolute cross section predic- 
tions. The measurements presented here include 
differential cross sections that may be useful con- 
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straints upon parton distribution functions. 

We use a coordinate system with the z 
axis along the proton beam, transverse coor- 
dinate perpendicular to the beam, azimuthal 
angle <f>, polar angle 8, and pseudorapidity 
j] = — lntan((9/2). The analysis uses the CDF 
calorimeters [9], which cover the pseudorapidity 
range \rj\ < 4.2. The calorimeters are constructed 
in a tower geometry and are segmented in depth 
into electromagnetic and hadronic components. 
The calorimeter towers are 0.1 unit wide in 77. 
The tower widths in <f> arc 15° in the central re- 
gion and 5° for \r}\ greater than approximately 
1.2. 

We begin by considering events from the 
data sample selected by the trigger requirement 
Y^E T > 175 GcV. We refer to this 175 GeV as 
E^l below. Event reconstruction uses a cone al- 
gorithm [4] described in more detail below. The 
transverse energy is defined as Et = E sin 8, 
where E is the scalar sum of energy deposited in 
the calorimeter within a particular cone and 8 is 
the angle between the beam direction in the lab- 
oratory frame and the cone axis. All calorimeter 
energy clusters [4] in the event with Et > 10 GeV 
are summed. The three leading jets in the labo- 
ratory frame are used as the basis of transforma- 
tion into the three-jet rest frame. In the three-jet 
rest frame, the incoming partons are, by conven- 
tion [10], labeled partons 1 and 2, and their mo- 
menta are designated p\ andp2, respectively. The 
highest energy jets in this frame have energies la- 
beled -E3, E4, and P5 and are ordered such that 
E 3 > E4 > £5. The outgoing partons associ- 
ated with these jets are correspondingly labeled 
partons 3, 4, and 5. 

A three-jet system in the massless parton ap- 
proximation can be uniquely described by five in- 
dependent variables (see Figure 5 in [11]). We use 
the following: 

1. the invariant mass of the three-jet system, 
TO3.7 



2. the cosine of the angle #3 between the av- 
erage beam direction (pav = Pi — P2) and 



parton 3 in the three-jet rest frame: 
Pay ■ P3 



cos 8% 



IpavI 1^3 



3. the cosine of the angle ip* between the plane 
containing the average beam direction and 
parton 3 and the plane containing partons 
3, 4, and 5 in their center of mass frame: 

,* (P3 X PAV) ■ (P4 X p 5 ) 
\P3 X PAV I \Pi X p 5 \ 

4. the Dalitz variable A3 (see below) for the 
leading jet, and 

5. the Dalitz variable A4 (see below) for the 
next-to-leading jet. 

The invariant m 3 j is calculated by sorting jets 
by their energies in the laboratory frame, boost- 
ing to the rest frame of those with the three 
highest energies, re-sorting jets by energy in that 
frame, then computing 7713,7 = X^=3 E\, where 
the Ei are the energies of jets 3, 4, and 5 in the 
rest frame. We have investigated the probability 
that a jet with energy less than the weakest of 
the three jets in the laboratory frame may have 
an energy greater than E 5 in the 3-jet rest frame 
from which it is excluded by this algorithm. The 
restriction imposed by the cut on full trigger effi- 
ciency (see below) makes this probability negligi- 
ble. 

The Dalitz variables, A;, are defined as A; = 
2 • Ei/msj, (i = 3, 4, 5). Momentum conservation 
restricts the ranges of the Dalitz variables to 

§ < X 3 < 1, 

- < Xa < 1, and 
2 ~ 



< A 5 < 



A set of trigger and offline requirements [12] 
rejects events associated with cosmic rays, beam 
halo, and calorimeter malfunctions. Events 
are required to have a reconstructed primary 
vertex, defined as the vertex with the largest 
^ P; (where P is the total momentum of par- 
ticle i leaving the vertex in the event), within 
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\z\ < 60 cm. Events are denned to have resolved 
multiple interactions if a second vertex with at 
least ten associated tracks is reconstructed in the 
vertex track detector, and if that vertex is sep- 
arated from the primary one by at least 10 cm. 
Because multiple interactions can change the jet 
multiplicity in an event, for example, misiden- 
tifying two-jet events as three-jet events, events 
with resolved multiple interactions are removed. 
The number of events with unresolved multiple 
interactions in which an additional jet could be 
misidcntified is estimated to be less than 2% so 
no correction for them is applied. The result- 
ing effective total integrated luminosity of the 
data sample is 77±4 pb _1 , where the uncer- 
tainty reflects both the overall luminosity uncer- 
tainty (4.2%) and the uncertainty (0.5%) associ- 
ated with the removal of resolved multiple inter- 
actions. 

An iterative cone algorithm [4] with cone radius 
R = ^(Ary) 2 + (A</>) 2 = 0.7 is used to identify 
jets. Here Ar] = t]2 — f]i and A0 = fa — <f>i- The 
subscripts 1 and 2 correspond to the axes of the 
cone and calorimeter tower, respectively. 

Jets that share towers are combined if the to- 
tal Et of the shared towers is greater than 75% 
of the Et of either jet; otherwise the towers are 
assigned to the nearest jet. Jet energies are cor- 
rected [4] for errors in the absolute and relative 
energy scales and for additional energy associated 
with the underlying event. Since partons that are 
radiated out of the cone lead to the same losses in 
the theoretical calculation and in the data, out- 
of-cone corrections are not applied. The Et of a 
jet is calculated from the reconstructed position 
of the primary event vertex. All three leading jets 
are required to have Et > 20 GeV and \t]\ < 2.0. 
Events with fewer than three jets are rejected. To 
avoid collinear soft gluon instability in the itera- 
tive jet clustering algorithm [13], a cone overlap 
cut is imposed: events are rejected if the distance 
AR in r]-(f> space between the axes of any two of 
the three leading jets is less than 1.0 (see Figure 
5 of [4], which shows that this selection require- 
ment reduces to approximately zero the proba- 
bility of the two jets being merged by the clus- 
tering algorithm). To exclude regions in which 
the geometrical acceptance [11] is less than about 



95%, we require |cos^| < Jl~ (E$£ /c 2 /m 2 j) 2 , 
where mij is the mass of the two leading jets in 
the three-jet system and is defined analogously to 

We require full trigger efficiency, which occurs 
when ^2 3 j ots Et > 320 GeV, where the sum is 
over the three highest energy jets in the event 
with corrected E T > 20 GeV [14]. The data are 
compared to the theoretical prediction by sort- 
ing events into bins of size 0.02 x 0.02 in X3-X4 
space, the Dalitz plane. Figure 1 shows the Dalitz 
distribution of data that remain after all of the 
selection requirements have been applied. 
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Figure 1. The three-jet data, after all selection 
requirements have been applied. The energy cor- 
rection procedure (see text) has not been applied. 
The figures at the corners of the distribution rep- 
resent typical three-jet topologies in those regions 
of the Dalitz plane. 



Before the final binning is done, the data are 
corrected for the effects of the combination of de- 
tector resolution and energy mismeasurement. A 
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correction factor is determined for each bin in the 
plane as follows. A sample of events is gener- 
ated at the parton level with the HERWIG Monte 
Carlo. The final state partons are hadronized. 
The events are then binned in the Dalitz plane. 
The same events are next passed through the 
CDF detector simulation and rcbinned. For each 
bin the ratio of the number of events after and 
before detector simulation is computed. This ra- 
tio (ranging from 0.85 to 1.5) is the factor subse- 
quently used to correct the number of events in 
each data bin. The data are also corrected for the 
z-vertex cut efficiency and then normalized to the 
effective total luminosity. 

The principal sources of systematic uncer- 
tainty [15] on the cross section are those on 
the absolute and relative (^-dependent) jet en- 
ergy scales. The uncertainty on absolute jet en- 
ergy derives from the resolution on the calibra- 
tion of the calorimeter (uncertainty 1.3%-1.8%, 
and i?T-dependcnt), the uncertainty associated 
with choice of jet fragmentation model (decreas- 
ing from 1.7% to 1.2% with increasing &r), the 
uncertainty associated with calorimeter stability 
over time (1%), and the uncertainty on the correc- 
tion for the contribution of the underlying event 
(1 GeV). The uncertainty on the relative jet en- 
ergy scale ranges from 2% to 6%. Uncertainties 
are also associated with the measurement of the 
effective total integrated luminosity (4.2%) and 
with the z-vcrtex cut efficiency (2%). There is 
also an uncertainty of less than 5% associated 
with the implementation of simulated events in 
the correction procedure. The upper (lower) lim- 
its on these uncertainties are added (subtracted) 
from the four-momenta of the jets in the data 
sample to obtain the systematic uncertainties on 
the cross section associated with each contribu- 
tion. The uncertainties are then combined to 
produce the total experimental systematic uncer- 
tainty. 

The Trirad calculation consists of 2 — > 3 parton 
processes at one loop and 2-^4 parton processes 
at tree level. For gluons g, incoming quarks q, 
and outgoing quarks Q or Q' , the subproccsscs 
involved are gg -» ggg, qq -> ggg, qq -> QQg, 
and those related by crossing symmetry, all com- 
puted to one loop; and gg -> gggg, qq -> gggg, 



qq — > QQgg, and qq — > QQQ'Q' and the crossed 
processes computed at tree level. The program 
uses the "subtraction improved" [13] phase space 
slicing method to implement infrared cancella- 
tion. 

The cross section is predicted with the 
CTEQ4M [16] parton distribution function 
(PDF) for each bin in the Dalitz plane. The result 
is multiplied by the effective total integrated lu- 
minosity of the data to predict a number of events 
in each bin. We restrict the prediction to bins for 
which X% < 0.98; this is necessary as the pertur- 
bative expansion is not reliable where the three- 
jet configuration approaches a two-jet configura- 
tion. The comparison between the data and the 
calculation is made for 215 bins. 

Figure 2 compares Dalitz distributions of the 
data and the absolute theoretical prediction. The 
theoretical distribution is more strongly peaked — 
a trend that persists in comparisons with all 
members of the CTEQ4 family 2 of PDFs. This 
trend, in which the edges of the Dalitz plane are 
more populated by data than by the prediction, 
may give some indication of the size of the higher 
order contributions to the cross section. 

The data and theory are compared in two dif- 
ferent ways. In Figure 3, we compare the shapes 
of their Dalitz distributions by normalizing the 
data and theory predictions to the same num- 
ber of events. In Figure 4, we normalize theory 
to the experimental luminosity and compare the 
absolute values of the cross sections that are ob- 
served and predicted. In both figures, the predic- 
tion is made using the CTEQ4M parton distri- 
bution function, and the difference between ob- 
served and predicted number of events, scaled by 
the number of predicted events, is computed. 

The theoretical prediction for the cross sec- 
tion, using CTEQ4M and all bins in the Dalitz 
plane but those with X 3 > 0.98, is 473 ± 
2(stat.)±i| (scale) jJ£(PDF) pb. The theoretical 
uncertainty associated with choice of renormal- 
ization and factorization scales, /iR and /zp re- 



2 The CTEQ4 family includes CTEQ4A1, CTEQ4A2, 
CTEQ4M, CTEQ4A5, and CTEQ4A6, which differ in the 
value of a B input to their global fit, and CTEQ4HJ, for 
which a higher statistical emphasis was given to the high 
E T data from CDF. 
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Figure 2. The event density in the Dalitz plane 
for (a) the data and for (b) the prediction by 
the NLO Monte Carlo calculation with CTEQ4M, 
normalized to luminosity. 
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spectively, is estimated by varying the scales, 
whose default value is Et, to values of Et/2 and 
2i?x while keeping /iR = /ip. The theoretical un- 
certainty associated with choice of PDF is esti- 
mated from the spread in the predictions gener- 
ated with all members of the CTEQ4 family. The 
measurement is not sensitive to the value of a s as 
is also shown in [17]. The measured cross section, 
using all bins in the Dalitz plane but those with 
X 3 > 0.98, is 458±3(stat.)l^ 3 (syst.) pb. This 
is consistent with the theoretical prediction and 
with a previous CDF measurement [14] after cor- 
rections are made for the efficiencies of additional 
cuts introduced in this analysis. The measured 
cross section, using all bins in the Dalitz plane, is 
466±3(stat.)±7o 7 (syst.) pb. 

The measurements at high X 3 may provide use- 
ful constraints on future theoretical models in 
that region. It appears that up to NLO the theory 
predicts more soft radiation than the data have 
in the region where the primary partons are ap- 
proximately back-to-back. The data, especially 
in the region above X 3 = 0.98 where a perturba- 
tive expansion is not reliable, may be useful input 
to theoretical models of gluon-emission processes. 
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Figure 3. The fractional difference between 
the data and the theoretical prediction, using 
CTEQ4M, as a function of X4, averaged over X 3 . 
The vertical bands show the systematic uncer- 
tainties. The error bars show the statistical un- 
certainties in cases where those are larger than 
the size of the symbol used. The prediction is 
normalized to the data to facilitate a comparison 
of the shapes of the distributions. 
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Figure 4. The fractional difference between cor- 
rected data and the NLO prediction, using the 
CTEQ4M parton distribution function, as a func- 
tion of A3 for various X4 bins. Error bars re- 
flect statistical uncertainty for cases in which it is 
larger than the size of the symbol used. Shaded 
bands indicate systematic uncertainty. The pre- 
diction is normalized to the luminosity of the 
data. 
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